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High-repetition-rate pulses have widespread applications in the fields of fiber communications, frequency 
comb, and optical sensing. Here, we have demonstrated high-repetition-rate ultrashort pulses in an all-fiber 
laser by exploiting an intracavity Mach-Zehnder interferometer (MZI) as a comb filter. The repetition rate of 
the laser can be tuned flexibly from about 7 to 1 100 GHz by controlling the optical path difference between 
the two arms of the MZI. The pulse duration can be reduced continuously from about 10.1 to 0.55 ps with 
the spectral width tunable from about 0.35 to 5.7 nm by manipulating the intracavity polarization 
controller. Numerical simulations well confirm the experimental observations and show that filter-driven 
four- wave mixing effect, induced by the MZI, is the main mechanism that governs the formation of the 
high-repetition-rate pulses. This all-fiber-based laser is a simple and low-cost source for various applications 
where high-repetition-rate pulses are necessary. 

Ultrafast lasers with repetition rates in the megahertz and gigahertz regimes and beyond have been inves- 
tigated for their immense applications in the fields of fiber optic communications 1 " 4 , materials proces- 
sing 5,6 , frequency comb 7 , and nonlinear optics 8 " 13 . In particular, ultrafast fiber lasers exhibit inherent 
advantages such as high gain coefficient, wide operating wavelength range, excellent heat dissipation, and robust 
mode confinement 14 " 20 , and offer an attractive platform for implementing high-repetition-rate (HRR) pulse 
sources. In the early demonstrations of achieving HRR pulse sources, an active mode-locking scheme was 
typically proposed, e.g., a fiber laser with repetition rate of 10 GHz was reported 21 . However, further enhancement 
of the repetition rate is limited by the available electrical bandwidth of the modulator and its driving electron- 
ics 21,22 . Currently, various approaches have been exploited to generate HRR pulses in passively mode-locked fiber 
lasers, which are simpler in construction and are free from electrical bandwidth limitation. Therefore, pulses with 
higher repetition rate could be achieved in passive mode-locking, compared with active mode-locking. For 
example: in fundamentally mode-locked fiber lasers, pulse repetition rates in the range 10-20 GHz have been 
obtained by minimizing cavity length and employing thin film absorbers based on nanotubes and graphene 23,24 . 
However, further scaling to higher repetition rate is limited by the physical length of the fiber-based cavity. For > 
1000 GHz, the necessary cavity length L should be less than 100 urn, assuming a typical refractive index n = 1.46 
for a typical single-mode-fiber based linear cavity. For such a short L, it is not easy to provide enough gain for 
lasing and also offering reasonable output power (e.g., > 1 mW). For harmonic mode-locking, where multiple 
pulses circulate simultaneously in the cavity per round- trip time 4,25,26 , repetition rate up to 10 GHz has been 
realized and such a laser can operate at as high as the 634 th harmonic of the fundamental cavity frequency 27 . But 
the harmonic mode-locking operation of laser is hardly controllable in practice 4,25 . 

Recently, HRR pulse trains (up to 200 GHz) were generated in a fiber laser incorporating an intracavity Mach- 
Zehnder interferometer (MZI) 28 . However, a semiconductor saturable absorber mirror (SESAM) was employed 
to achieve the mode-locking operation and the pulse duration appeared to be independent of the interferometer 
settings. It is worth noting that Yoshida et al. demonstrated a low- threshold 1 15-GHz pulse train in a long-cavity 
fiber laser without any saturable absorber 29 . This scheme was later reinterpreted in terms of dissipative four- wave 
mixing (DFWM) where only two frequencies experience positive gain and transmit their energy by four-wave 
mixing (FWM) to high-order harmonics that undergo linear loss owing to the bandpass nature of the filter 30,31 . 
Peccianti et al. have recently realized stable HRR pulses in a mode-locked fiber laser based on a high-Q factor 
microring resonator 32 . They term this new mode-locking scheme filter-driven four-wave mixing (F-D FWM) 
because the microring resonator works as both a nonlinear element and a Fabry-Perot filter 32,33 . Despite the 
elegance of the FWM-based methods, the repetition rates of lasers based on these schemes are typically 
fixed 29 " 32,34 , and even researchers have to overcome the difficulty of splicing the fiber and microcavity 32,33 . 
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Figure 1 | Setup of the HRR fiber laser. LD, laser diode; WDM, 
wavelength division multiplexer; EDF, erbium-doped fiber; SMF, single- 
mode fiber; PC, polarization controller; OC, optical coupler; OTDL, 
optical time delay line. The MZI is composed of two 5 : 5 OCs and an 
OTDL. 



Therefore, it is of great significance to study and realize repetition- 
rate-tunable duration- controllable ultrafast fiber lasers with an all- 
fiber approach. 

Here, we propose a compact all-fiber laser with HRR pulses, in 
which an intracavity MZI plays a crucial role in the realization of the 
mode locking. By adjusting the optical path difference between the 
two arms of the MZI, the repetition rate of the pulses can be flexibly 
tuned from ~ 7 to 1 100 GHz. Meanwhile, the spacing between adja- 
cent peaks of the spectrum increases from ~ 0.06 to 9.5 nm. 
Interestingly, the pulse duration can be adjusted continuously from 
~ 10.1 to 0.55 ps. The simulation results closely reproduce experi- 
mental observations. Our study provides a simple and low- cost 
method for generating HRR pulses by utilizing an all-fiber MZI. 

Results 

Setup of the HRR fiber laser. The schematic diagram of the laser 
setup is shown in Fig. 1. The resonator consists of a 17.6 m erbium- 
doped fiber (EDF: CoreActive C300) with a 3 dB/m absorption at 
980 nm, an isolator, two polarization controllers (PCs), and a MZI. 
The other fibers are single-mode fibers (SMFs) with a total length of 
9.4 m. The dispersion parameters D for EDF and SMF are about 
— 9 ps/(nnrkm) and 17 ps/(nnrkm) at 1550 nm, respectively. 
The net dispersion of laser cavity is about —0.0018 ps 2 , which 
approaches zero, and then the Gaussian-profile pulse will be 
generated from the laser system. The MZI is constructed by 
splicing two identical 5 : 5 optical couplers (OCs). A tunable optical 
time delay line (OTDL) with a pigtail of 1.5 m is placed in one of the 
arms and the total length between two OCs is ~ 2 m. The MZI works 
as a spectral filter for wavelength selection, and the transmission 
spectrum can be precisely controlled by the OTDL. The output 
ratio of the OC is 50%. Two 980-nm laser diodes (LDs) and 
wavelength division multiplexers (WDMs) are used to provide 
bidirectional pumps. The total pump power is tunable from 0 to 
1 W, which is sufficient to induce strong nonlinear effects in 
cavity. The isolator is used to ensure unidirectional operation in 
the oscillator. Note that HRR pulses can be generated by using 
either a polarization- sensitive isolator or a polarization-insensitive 
isolator in the experiment and output pulses are almost independent 
of the two different isolator settings. 

Experimental observations. Self- starting HRR pulse trains can 
always be obtained in the fiber laser by precisely adjusting the set- 
tings of PCs and OTDL. When the forward and backward pumps are 
increased to 400 mW, the average output power is ~ 50 mW and is 
insensitive to the time delay of OTDL. Figures 2(a)-(d) show the 
output spectra for optical path differences of 12.27, 6.91, 3.4, and 
0.27 mm, respectively. The optical spectra exhibit comb profiles with 
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Figure 2 | (a)-(d) Optical spectra of the pulses and (e)-(h) corresponding 
AC traces at optical path differences AL of 12.27, 6.91, 3.4, and 0.27 mm, 
respectively. 

spectral separations of 0.193, 0.37, 0.75, and 9.5 nm, and the 
bandwidths of the spectral envelopes are estimated as 1.66, 1.6, 
1.62, and 9.5 nm, respectively. The slight asymmetry of the spectra 
may be attributed to the uneven gain characteristics of the EDF. The 
relevant autocorrelation (AC) traces, as shown in Figs. 2(e)-(h), 
demonstrate that the laser produces equidistant pulse trains with 
temporal spacings of 41, 23, 11, and 0.9 ps, corresponding to pulse 
repetition rates of about 24, 43, 91, and 1100 GHz, respectively. If a 
Gaussian profile is assumed for fitting, the durations of the individual 
pulses in Figs. 2(e)-(h) are estimated as 2.20, 2.23, 2.21, and 0.37 ps. 
Then, the time bandwidth products (TBPs) of the pulses are 
calculated as 0.45, 0.44, 0.44, and 0.43, indicating that the output 
pulses are nearly chirp-free (i.e., transform-limited). Practically, 
the highest repetition rate of the pulses is limited by the available 
gain bandwidth of EDF. For example, only two spectrum peaks are 
generated within the gain bandwidth, as shown in Fig. 2(d). The HRR 
pulses with well- separated spectral peaks can be used for wavelength- 
division-multiplexed applications, an approach that may be 
economically advantageous over employing a large number of 
individually selected continuous-wave sources. 

Figure 3 shows the spectral separation as a function of the optical 
path difference AL. Based on the transmission function of MZI, the 
spectral separation Ak of adjacent peaks is given by Ak = k 2 lAL (see 
Methods for details), where X is the operating wavelength. The repe- 
tition rate of the pulses is equal to the frequency difference of adja- 
cent spectral peaks and is expressed as / = cAX/X 2 , where c is the 
speed of lightwave in vacuum. Experimental results show that the 
spectral separation increases from ~ 0.06 to ~ 9.5 nm by decreasing 
AL. Correspondingly, the repetition rate / of the pulses increases 
from ~ 7 to ~ 1100 GHz. It is evident that the wavelength spacing 
and repetition rate of the output pulses can be flexibly tuned over a 
large range with the OTDL in the MZI, which is a distinct advantage 
over the fiber laser mode locked with the micro-resonator 32 . 

A notable feature of the laser is that the spectral width and pulse 
duration could be varied by adjusting the settings of the intracavity 
PCs if a polarization- sensitive isolator is used. Here, the optical path 
difference of the interferometer is fixed at 12.27 mm, and thus the 
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Figure 3 | Spectral separation versus optical path difference AL. The 

repetition rate f=c Ak/k 2 . The theoretical curve of spectral separation is 
obtained from the equation Ak = k 2 /AL. 

spectral separation and pulse-pulse separation keep unchanged at 
0.193 nm and 41 ps respectively in the following experiment. As 
shown in Figs. 4(a)-(h), the spectral widths of pulses are about 
0.35, 0.74, 0.91, and 5.7 nm, and the corresponding pulse durations 
are about 10.1, 4.9, 3.7, and 0.55 ps, respectively. Although the pulses 
have different spectral widths and durations, they almost process the 
same TBP of about 0.44. The small periodic modulation on the 
spectral pedestal can be attributed to spectral filtering effect induced 
by the MZI. One can see from Fig. 4 that a broader optical spectrum 
corresponds to a narrower pulse, which coincides with the theory of 
soliton mode locking 35 . Furthermore, we find that the central wave- 
length of the spectrum can be adjusted from 1530 to 1570 nm by 
changing the polarization state. Here, the tuning range of the laser is 
practically limited by the available pump power and gain bandwidth 
of the EDF. The tuning feature of the pulse is possibly caused by the 
birefringence-induced spectral filtering effect. The PC-induced 
birefringence may introduce a wavelength- dependent polarization 
that, in conjunction with the polarization -sensitive isolator, will 
impose amplitude modulation on the optical spectrum. As a result, 
the transfer function of the cavity relies on the PC states, and mode 
locking tends to be realized at the wavelength with a higher trans- 
mission. This assumption can be further confirmed from Figs. 4(d) 
and 4(h), in which the optical spectrum exhibits modulation and the 
AC trace has small satellites. It is worth noting that the output spec- 
trum in Fig. 4(d) is distinct from that of DFWM mode locking 29 " 31 
where the spectral intensity decreases from center to edges signifi- 
cantly. The results indicate that the dissipative process plays a less 
important role in the HRR pulse formation. 

Numerical results. Figure 5 shows the numerical results of the HRR 
fiber laser. In the simulation, the initial signal is a weak light field. 
Due to the combined effects of spectral filtering, FWM, cavity 
dispersion, and fiber nonlinearity, the optical field evolves to a 
steady state and becomes symmetric after ~ 600 round trips. As 
shown in Fig. 5(a), multi- wavelength spectrum is gradually formed 
due to spectral filtering effect, and the central wavelengths transfer 
their energy to side wavelengths through FWM effect. Figure 5(b) 
shows that the light field is splitted into several sub -pulses and 
eventually evolves to a stable state with the increase of round-trip 
number. One can see that the output pulses are uniform and are 
equally spaced with a separation of 1 1 ps in the temporal domain. 

The red curve in Fig. 6(a) shows the transmission spectrum of the 
MZI, which is characterized by a series of equally spaced transmis- 
sion peaks in the frequency domain. As shown in the inset of 
Fig. 6(a), the full width at half maximum (FWHM) of a single trans- 
mission peak is about 0.35 nm. The output spectrum at round trip of 
1200 (the blue curve in Fig. 6(a)) exhibits a multi-channel structure 
with a wavelength spacing of 0.75 nm. We observe that lasing can 
only take place at the transmission peaks of the MZI, and the 
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Figure 4 | (a)-(d) Optical spectra of pulses and (e)-(f) corresponding AC 
traces at different PC states. The optical path difference of the 
interferometer is fixed at 12.27 mm, and the spectral separation keeps 
unchanged at 0.193 nm. If a Gaussian profile is assumed for fitting, the 
pulse duration is the width of AC trace divided by a factor of ~ 1.41. 



wavelength spacing of the comb-shape spectrum is equal to that of 
the transmission spectrum. The corresponding temporal profile and 
AC trace of the HRR pulses are shown in Figs. 6(b) and 6(c), respect- 
ively. The calculated AC trace is similar to that of the experimental 
results in which the intensity of AC peaks decreases slightly from 
center to edges. Our numerical results have fully confirmed experi- 
mental observations in Figs. 2(c) and 2(g), and demonstrated that the 
MZI -induced filtering effect dominates the pulse formation. 

Discussion 

In the experiments, the measured FWHM A/I FWH m of a single spec- 
tral peak ranges from 0.015 to 0.025 nm, which is limited by the 
resolution of optical spectrum analyzer (0.015 nm). For example, 
the AA FWH m of a single peak in Fig. 2(c) is about 0.02 nm. The 
longitudinal mode spacing Av is ~ 7.61 MHz, which corresponds 
to spectral separation Ak LM of ~ 6.08 X 10" 5 nm. Then there are 
about 300 modes in a single spectral peak. In comparison with a 
traditional mode-locking fiber laser that contains ~ 4 X 10 5 modes 
(A2 FWHM « 0.96 nm, Av « 0.29 MHz, Ak LM « 2.3 X 10" 6 nm) 16 , 
the modes here are much less in a single spectral peak. As a result, the 
fewer modes (—300) in this experiment are more easily synchronized 
than the more modes (~4 X 10 5 ) in Ref. 16 from a theoretical point 
of view 35 . 

Because every spectral peak is composed of many longitudinal 
modes, the intensity of the output pulses can be modulated with 
~ 7.61 MHz, corresponding to ~ 131 ns. Due to the limited oper- 
ating range of the autocorrelator (100 ps), the modulation period can 
not be directly obtained from the AC trace. We have theoretically 
investigated the modulation of pulse trains by summing the longit- 
udinal modes in every spectral peak 35 . In the calculation, the mode 
separation and the optical path difference are assumed as 7.61 MHz 
and 3.4 mm, respectively, which are coincided with that of the 
experiment in Fig. 2(c). The phase of each mode is assumed to change 
randomly from 0 to 1. The theoretical results show that the HRR 
pulses are modulated with a period of ~ 131 ns. The calculated AC 
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envelope decreases slightly from left to right, which is similar to that 
of Fig. 2(g). 

We have replaced the polarization- sensitive isolator by a polariza- 
tion-insensitive isolator in the cavity. HRR pulses with multi- wave- 
length spectrum are still achieved, and the repetition rate can be 
tuned over a wide range by the OTDL. The output spectrum and 
AC trace are similar to that of Fig. 2. However, the central wavelength 
and bandwidth of the spectrum are hard to be adjusted by PCs as 
other cavity components are polarization-insensitive. Our above 
experimental observations (i.e., successful demonstration of mode- 
locking using a polarization-insensitive or polarization-sensitive 
isolator) immediately demonstrate that the mode-locking is not 
achieved by nonlinear polarization rotation technique. This is further 
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Figure 6 | (a) Transmission of the MZI and output spectrum of pulses, (b) 
pulse train, and (c) AC trace at round trip of 1200. The inset in Fig. 6(a) is 
the transmission spectrum in linear scale. 



confirmed by that removing MZI from cavity cannot lead to mode- 
locking with the polarization -insensitive isolator. Compared with the 
laser in Ref. 28 that a SESAM is inserted in the cavity as a mode- 
locking element, only a MZI is used to generate HRR pulses in this 
scheme. Our results confirm that the F-D FWM induced by the MZI 
mainly contributes to the generation of pulses. 

The formation of the HRR pulses may be understood as follows. 
During the pulse amplification process in EDF, the optical spectrum 
broadens strongly with the increase of pulse intensity due to the SPM 
and other nonlinear effects. With the help of the intracavity MZI, the 
optical spectrum of the pulses is modulated periodically. And then, 
the energy in the cavity is confined in several selective wavelengths. 
Due to the FWM and other nonlinear effects, the high-intensity 
central wavelengths distribute their energy to side wavelengths and 
form the comb-shape spectrum. Because each wavelength has the 
regular phase and position during the FWM process, they interfere 
and form the HRR pulses in the temporal domain. This process 
occurs at each circulation, and stable HRR pulses are formed when 
the optical field becomes self-consistent. Compared with methods in 
Refs. 29-32, only a fiber-based MZI is used as a spectral filter to 
achieve HRR pulses in our cavity. As the spectral filtering effect plays 
a key role in the FWM process, we attribute the pulse formation to F- 
D FWM effect. 

Methods 

Measurement method. An optical spectrum analyzer (Yokogawa AQ-6370) and a 
commercial autocorrelator (APE PulseCheck) are employed to characterize the laser 
output simultaneously. The repetition rate is deduced from the pulse-pulse separation 
measured by the autocorrelator. 

Numerical simulation. To get an in-depth understanding of the experimental 
observations, we numerically simulate the dynamic evolution of HRR pulses based 
on the nonlinear Schrddinger equation (NLSE) 35 ' 36 . We follow the propagation of the 
optical pulses in the laser cavity and consider functions of every cavity components 
on pulses. When pulses meet a cavity component, we then discretely multiply the 
light field by the relevant transformation matrix. After one round-trip circulation in 
the cavity, we use the obtained results as the input of the next round of calculation 
until the optical field becomes self-consistent. As the pulse duration is much 
larger than the delayed Raman response, the term related to Raman interaction 
can be neglected 35 " 38 . Pulse propagation within each fiber in laser cavity is modeled 
by 36 ' 38 



du(z,t) k 2 d 2 u(z,t) 

dz + *T dt 2 



-iy\u(z,t)\ u(z,t) 



g S 2 u(z,t) 
2Ql dt 2 ' 



(1) 



where u(z,t) designates the slowly varying envelope of the optical field. The variables 
t and z represent the time and propagation distance, respectively. k 2 is the second- 
order dispersion coefficient, and y is the cubic refractive nonlinearity of the fiber. Q g 
is the gain bandwidth of the laser, g is the saturable gain of the fiber, which is set as 
zero for SMF. For the EDF, the gain saturation is expressed asg = g Q 'exp( — E p /E s ) 39 , 
where g 0 , E p and E s are the small- signal gain coefficient, pulse energy, and gain 
saturation energy, respectively. There is not any saturable absorber in the numerical 
model. 

The filtering effect induced by the MZI is taken into account in our simulations. 
The transfer function of the interferometer can be expressed as 35 : 
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T=[l-cos(* 1 + * M .)]/2. 



(2) 



And the linear and nonlinear parts of the relative phase shift are given by 

® L = P 0 AL = 2nAL/A, O nl = yP 0 AL/2, (3) 

where /3 0 is propagation constant, y is the nonlinear coefficient, P 0 is the pulse peak 
power, AL is the optical path difference between the two arms of the MZI, and X is the 
lasing wavelength. As O nl << O l in this experiment, the simplified transfer function 
of interferometer can be described as 



T=[l -cos(2ttAL/4)]/2. 



(4) 



Since the linear phase shift is frequency-dependent, the transmission spectrum of the 
interferometer depends on the wavelength of light and the modulation frequency is 
mainly determined by AL. The modulation period corresponds to <J> L = 2%. The 
relation between modulation period and spectral separation AX is deduced as follows: 



2nAL/A 1 — 2%AL/X 2 = 2%, 



AXAL/(X l X 2 ) = l. 



(5) 



(6) 



Here, AX = X 2 — X\. Because AX is very small, X\X 2 approximately equals X 2 . Then, the 
spectral separation between adjacent transmission peaks can be expressed as: 



AX = X 2 /AL. 



(7) 



The NLSE is solved with a predictor-corrector split- step Fourier technique 38,40 . The 
gain spectrum of EDF is assumed to have a Gaussian profile with a 3-dB bandwidth of 
18 nm. We use the following parameters to match the experimental conditions, c = 3 
X 10 8 m/s,£ 0 = 3 dB/m, n = 1.46, AL = 3.4 mm, a = 0.2 dB/km, Q g = 18 nm, E s = 
15 pj.y = 3 W^-km^andA^ = 11.47 ps 2 /km for EDF; y = 1 W^-km" 1 and k 2 = 
-21.67 ps 2 /km for SMF. The lengths of EDF and SMF are 17.6 and 9.4 m, respect- 
ively. The other parameters are the same as that of experiments. 
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